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Abstract 

Theories that spontaneously break Lorentz invariance also violate diffeo- 
morphism symmetries, implying the existence of extra degrees of freedom 
and modifications of gravity. In the minimal model ( "ghost condensation" ) 
with only a single extra degree of freedom at low energies, the scale of 
Lorentz violation cannot be larger than about M ~ 100 GeV due to an 
infrared instability in the gravity sector. We show that Lorentz symmetry 
can be broken at much higher scales in a non-minimal theory with addi- 
tional degrees of freedom, in particular if Lorentz symmetry is broken by 
the vacuum expectation value of a vector field. This theory can be con- 
structed by gauging ghost condensation, giving a systematic effective field 
theory description that allows us to estimate the size of all physical ef- 
fects. We show that nonlinear effects become important for gravitational 
fields with strength > g, where g is the gauge coupling, and we ar- 
gue that the nonlinear dynamics is free from singularities. We then ana- 
lyze the phenomenology of the model, including nonlinear dynamics and 
velocity-dependent effects. The strongest bounds on the gravitational sec- 
tor come from either black hole accretion or direction-dependent gravita- 
tional forces, and imply that the scale of spontaneous Lorentz breaking is 
M < Min(10^^ GeV, g"^ 10^^ GeV). If the Lorentz breaking sector couples 
directly to matter, there is a spin- dependent inverse-square law force, which 
has a different angular dependence from the force mediated by the ghost 
condensate, providing a distinctive signature for this class of models. 



1 Introduction 



As the 100th anniversary of special relativity comes to a close, Lorentz invariance 
reigns as arguably the most important symmetry in modern physics. A good way 
to understand the central role of Lorentz invariance is to consider how it could be 
violated. Like baryon and lepton number in the standard model, Lorentz invariance 
could be an accidental symmetry of the leading interactions in quantum field theory. 
Alternatively, Lorentz symmetry could be one limit of a more fundamental symmetry, 
just as the Galilean group is the small velocity limit of the Lorentz group. But in the 
absence of experimental data, it is difficult to guess what deeper organizing principle 
could replace Lorentz invariance. Indeed, the Lorentz group itself is one of many 
deformations of the Galilean group, and it was the experimental crisis of the ether 
that drove the transition from absolute time to the space-time continuum in 1905. 

Here, we consider the less radical possibility that Lorentz invariance is a good 
symmetry at high energies, but is spontaneously broken at an energy scale M. Some 
of the consequences of breaking Lorentz invariance have been extensively explored 
in the literature. If the symmetry breaking sector couples to standard model fields, 
there will be Lorentz-violating operators in the effective theory that give rise to CPT 
violation and a variety of preferred frame effects (see Ref. [1, 2, 3, 4, 5, 6, 7] and 
references therein). A second aspect that has been received less attention is that 
spontaneous breaking of Lorentz invariance implies the existence of new degrees of 
freedom analogous to the Goldstone bosons that arise from spontaneous breaking of 
internal symmetries. 

The minimal model of this kind was analyzed in Ref. [8] , which contains a single 
Goldstone degree of freedom. It can be viewed as the result of ghost condensation 
(in the same way that ordinary spontaneous symmetry breaking can be viewed as 
"tachyon condensation" ) and so we call the extra degree of freedom a "ghostonc" bo- 
son. If the ghostone mode couples to the standard model, then it gives rise to exotic 
spin- and velocity-dependent forces [8, 9]. In a gravitational theory, breaking Lorentz 
invariance immediately implies a violation of diffeomorphisms, the gauge group of 
Einstein gravity, and thus the ghostone mode will mix with the metric. This is the 
Higgs mechanism for gravity, and is analogous to the familiar Higgs mechanism for 
gauge theories, where the mixing between a gauge field and a Goldstone boson gives 
rise to a massive spin-1 boson. Refs. [8, 10] analyzed the modification of gravity that 
arises in this way, and found the limit M < 100 GeV because of an infrared gravita- 
tional instability analogous to the Jeans instability for ordinary matter. Therefore, 
spontaneous Lorentz breaking is directly tied to modifications of gravity, and we can 
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place strong constraints on Lorcntz breaking by looking at gravitational physics. 

In this paper, we consider a non-minimal model of spontaneously broken Lorentz 
invariance that has three extra degrees of freedom instead of a single ghostone mode. 
We will show that in this model it is possible to break Lorentz invariance sponta- 
neously at much higher scales, as high as M ~ Min(10^^ GeV, GeV). The 
model can be thought of as arising from a vector field order parameter, and has been 
previously considered in the literature [11, 12, 13, 14, 15, 16, 17, 18, 19, 20, 21, 22, 23]. 
There are several features of our analysis that are new. First, we show that this mo- 
del can be viewed as a gauging (in the ordinary sense) of ghost condensation. This 
considerably clarifies the nature of the effective theory, and allows us to consistently 
estimate the size of all operators in a low energy expansion. Second, we consider 
the couphngs of the Goldstone degrees of freedom to matter, and show that there 
are new spin- and velocity-dependent forces with a different signature from the min- 
imal model. Third, we show that the theory has important nonlinear corrections in 
strong gravitational fields. We find that the nonlinear effects are much less dramatic 
than in ghost condensation (see Ref. [10]), and argue that the theory is free from 
caustic singularities. Finally, we consider the modification of gravity in this model, 
studying observable consequences such as velocity-dependent effects (such as gravita- 
tional Cerenkov radiation), nonlinear dynamics, and the accretion of the condensate 
by black holes. Our conclusion is that the model is consistent with all experimental 
limits for M < Min(10-^^ GeV, g"^ 10^^ GeV), where g is the gauge coupfing and other 
dimensionless coefficients are taken as 0{1). 

This paper is organized as follows. In Sec. 2, we describe the effective theory 
by directly introducing the extra degrees of freedom required by diffeomorphism in- 
variance. This is the most direct route to the effective theory and the modification 

of gravity. In Sec. 3, we describe the same effective theory as the gauged version 
of ghost condensation. This formulation makes it manifest that this is a theory of 
spontaneous breaking of Lorentz symmetry, and also makes it easy to understand the 
general power counting and the construction of the nonlinear theory. In Sec. 4, we 
couple the gauged ghost condensate to gravity and perform a linear analysis. We find 
that the modification of gravity is mild in this case and it is the form of a direction- 
dependent Newtonian potential. In Sec. 5, we study the nonlinear effects. In contrast 
to the ungauged ghost condensate, the would-be caustics are smoothed out by the 
gauge interaction and there are no dangerous instabilities. In Sec. 6, we consider 
gauged ghost condensate surrounding a black hole and calculate the rate of increase 
of the black hole mass due to accretion of the condensate. In Sec. 7, we compute 
the energy loss due to the ether Cerenkov radiation for a source moving faster than 
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the ghostone mode. In Sec. 8, we study various Lorentz-violating effects and tlieir 
constraints if the gauged ghost condensate couples directly to the Standard Model 
fields. Sec. 9 contains our conclusions. Some more detailed analyses and discussion 
are included in the appendices. 



In a non-gravitational theory, Lorentz symmetry can be broken explicitly without 
introducing any Goldstone degrees of freedom. The reason is that Lorentz symme- 
try is a physical symmetry that arranges matter states into multiplcts, so violating 
Lorentz symmetry will only affect the allowed couplings between matter states.^ As 
we will see, however, breaking Lorentz invariance implies a breaking of diffeomor- 
phism invariance, the gauge symmetry of Einstein gravity. Since a gauge symmetry 
is a redundancy of description rather than a physical symmetry, breaking diffeomor- 
phisms necessarily introduces additional degrees of freedom. The key point is that 
these new modes nonlinearly realize diffeomorphisms (as well as nonlinearly realiz- 
ing Lorentz transformations), and therefore their couplings to gravity are uniquely 
determined by diffeomorphism invariance. 

In this section, we show how to construct the effective theory below the Lorentz 
breaking scale by explicitly adding these degrees of freedom to the Lagrangian. We 
will also see that different ways of nonlinearly realizing diffeomorphisms give rise 
to different low energy effective theories. We review the case of ghost condensation 
which has only one new degree of freedom and introduce "gauged ghost condensation" 
which has three new degrees of freedom. We find that the mixing between gravity 
and the gauged ghost condensate is mild, which explains why the scale of spontaneous 
Lorentz breaking can be much higher in this theory. Other interesting generalizations 
are also possible, see for example Refs. [24, 25, 26]. 

The simplest way to break Lorentz invariance is to include the time component 
of a vector field in the Lagrangian: 



The arguments below are easily generalized to components of higher tensors. Eq. (2.1) 
breaks Lorentz symmetry down to 5*0(3) rotations, and also breaks diffeomorphism 
symmetry. Moreover, there is no way to couple A>C to the metric g/^^, to restore 

"'^Of course, one can spontaneously break Lorentz symmetry in a non-gravitational theory, which 
would give rise to Goldstone modes. The point is that in the absence of gravity, Lorentz violation 
can be explicit, whereas in a gravitational theory, Lorentz violation must be spontaneous. 



2 Lorentz Violation and Diffeomorphisms 
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diffeomorphisms. To see this, we work in the hnearized theory obtained by expanding 
the metric about flat space (where Lorentz invariance has an unambiguous meaning) 

Qixi^^Vixi^ + Ki^- (2-2) 

Under the infinitesimal diffeomorphism generated by ^ — the metric and 
transform as 

= d^i^u + (2-3) 

= ^PdpJ^ - dp^JP. (2.4) 

We have 

5{V^J')--V^d,ej^ (2.5) 

which makes it clear that there is no way to introduce factors of h^i, to make A>C 
invariant under diffeomorphisms. 

However, we can formally restore diffeomorphism invariance by introducing a 
scalar field tt that transforms as 

57r = eo- (2.6) 

This is the degree of freedom in the "ghost condensate" model of Ref. [8], and tt 
naturally has dimensions of length. The combination 

A£ = v^(J° + 9^7rJ'^) (2.7) 

is fully diffeomorphism invariant.^ Note that diffeomorphism invariance is realized 
nonlinearly, since the vr field transforms inhomogeneously. 

Because vr is a real propagating degree of freedom, we need to introduce a kinetic 
term for tt, but diffeomorphism invariance requires the kinetic term to depend on h^i,. 
Specifically, the time kinetic term 

-^Trkinetic = +|M^(7r — |/loo)^ (2.8) 

is diffeomorphically invariant in the linearized theory. Note that the sign is fixed 
by the requirement that tt fluctuations have positive energy. Eq. (2.8) shows that 
diffeomorphism invariance requires mixing between the new degree of freedom and 
gravity. In a gauge where tt = 0, Eq. (2.8) becomes a "wrong sign" mass term 

^In fact, Eq. (2.7) is invariant in full nonlinear Einstein gravity. 
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for hoo, suggesting an oscillatory Newtonian potential at distances larger than the 
inverse "mass" Tc ~ Mpi/M^. A more careful analysis shows that the potential is 
indeed modified at distances larger than Tc, but only on a much larger time scale 
tc ~ Mp^/M^ [8]. The theory is unstable on time scales of order tc, and demanding 
that this is longer than the present age of the universe gives M < 10 MeV, though 
nonlinear effects shghtly relax this bound to M ^ 100 GeV [10]. 

The strong bound on M arises because the way we have restored diffeomorphism 
invariance allows a graviton "mass" To find a theory that allows larger values of M, 
we look for a way of restoring diffeomorphism invariance that does not necessitate a 
graviton mass term. To do this, note that 

5[v^(i + ihoo)A ^V^d^er, (2.9) 

where the sum over i does not include any additional minus signs. This shows that 
we can restore difformorphism invariance by introducing three fields Oj transforming 
as 

Sai = di^o- (2.10) 

Note that defined in this way is naturally dimensionless. Because these fields trans- 
form as the gradient of a diffeomorphism parameter, the time kinetic term involves 
mixing with derivatives of /i^^: 

>Cakinetic ~ +M^{ai - ■ (2.11) 

It is easy to see that all diffeomorphism-invariant combinations of and h^i, involve 
derivatives acting on the metric, so modifications of Newton's law are suppressed. We 
will analyze this theory in more detail starting in Sec. 3, and we will see that this 
construction allows much larger values for the scale M where Lorentz symmetry is 
broken. 

To see the connection between the above theory and gauging ghost condensation, 
note that the diffeomorphism transformation Eq. (2.10) is very similar to an ordinary 
U{1) gauge transformation 

= d^X. (2.12) 

In fact, we can promote Eq. (2.10) to an ordinary gauge transformation by introducing 
a field tt that transforms under gauge transformations and diffeomorphisms as 

5Ti = io- A. (2.13) 

^Strictly speaking, only the Newtonian potential gets a mass; the propagating spin-2 gravitons 
are still massless. 



5 



The combination 



S = + d^i: (2.14) 

is then gauge invariant and transforms under diffeomorphisms as 

5~a^ = 9^eo- (2.15) 

In "unitary gauge" vr = wc have hi = Cj, so the spatial components can be identified 
with the fields introduced in equation Eq. (2.10). The extra degree of freedom Oq is 
naturally massive, since nothing forbids the term 

>Caomass = +lM\ao - |/loo)' (2.16) 

which is a mass term for in unitary gauge. 

As we will explain in sections 3 and 4, Eq. (2.16) reduces to Eq. (2.8) when we 
take the gauge coupling to zero, which explains why the modification of gravity 
is suppressed in gauged ghost condensation. We find that for canonically normalized 
fields in unitary gauge, Eq. (2.16) becomes 

9 M 
Cmsss^^(gA'o-9c^'') , 3c = ^^^, (2.17) 

where is the canonically normalized gauge field, is the canonically normalized 
Newtonian potential, and g is the gauge coupling. When g is smaller than g^ 
the (negative) mass term dominantly affects the Newtonian potential, yielding large 
modifications to gravitational potentials as in ghost condensation. When g is larger 
than gc, the mass term acts mostly on the "Coulomb" potential A^, shielding Lorentz 
breaking from the gravitational sector and allowing the scale M in gauged ghost 
condensation to be raised as high as M < Min(10^^ GeV, g^ 10^^ GeV). 

A discussion of the decoupling limit of gauged ghost condensation which makes 
contact with earlier literature on Lorentz breaking appears in App. A. 

3 Gauged Ghost Condensation 

In this section, we will start from the ghost condensate and construct the effective 
theory by gauging a global symmetry. This will clarify the power counting, and also 
show how the theory reduces to the ghost condensate in a particular limit. 
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3.1 A Review of Ghost Condensation 



We first review the "ghost condensate" model of Ref. [8]. This is an effective theory 
of a real scalar field with a global shift symmetry 

0^0 + A. (3.1) 

The effective theory is assumed to depend on a single scale Af , which gives both the 
strength of self-interactions of the field and the cutoff of the effective theory. The 
novel feature of this sector is that it has a vacuum (in flat spacetime) 

{d,4>) = v„ (3.2) 

where is a constant timelike vector. This spontaneously breaks Lorentz invariance 
by establishing a preferred time direction. We can then choose a Lorentz frame where 
Vq = c, Vi = for i = 1,2, 3. In this frame we have 

(0) = ct. (3.3) 

It is convenient to take to have units of time, so that c is dimensionless and can be 
set to 1 by rescaling the time coordinate. 

Expanding about the vacuum 

0=(0)+7r, (3.4) 

we can construct a systematic effective field theory for the fluctuations tt in derivatives 
simply by noting that 

(3.5) 

The constant part of (9^0 should be kept to all orders since there is no momentum 
suppression. On the other hand, terms with more than one derivative acting on only 
give rise to derivatives acting on the fluctuations, so there is a well-defined expansion 
at low energies and only the lowest order terms are important. The most general 
quadratic effective Lagrangian with at most four derivatives acting on tt is 

>Cghost = M^Pix) + iM2(Qi(x)(n0)2 + g2(^)9^09>9^a,0n0 

(3.6) 

+ QsiX){dyd''<pd,d,<Pf) + ---, 
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where we have assumed a — > — reflection symmetry, and 

X = d''<pd^<l>. (3.7) 

Expanding this Lagrangian about the vacuum using Eq. (3.5), the leading La- 
grangian is 

>Cghost = iM^TT^ - iM^ (^q( VV)2 + /?7fVV - 7ir2) + O (7^^ (aV)^) . (3.8) 

Here we have rescaled the fields so that P"{1) — |, which fixes the coefficient of the 
TT^ term, and used the fact that P'{1) = in the vacuum so that there is no tadpole 
term. The coefficients of the other terms are given by 

a = -Qi(l), 

/?-2gi(i) + g2(i), (3.9) 

7 = Qi(i) + Q2(i) + g3(i). 

We can further rescale tt and M to simplify this Lagrangian, but we will not make 
use of this freedom for now. We can see that there is no normal spatial kinetic term 
(Vtt)^, which implies the special dispersion relation for tt 

= «^ (3-10) 

to leading order in cu and k. 

3.2 Gauging Ghost Gondensation 

We now couple the theory above to a gauge field by gauging the shift symmetry 
in Eq. (3.1). That is, we promote Eq. (3.1) to a local transformation and introduce 
a gauge field to make the Lagrangian invariant. It is convenient to define the gauge 
field as — Ma^ so that the gauge field has dimensions of mass. We then write the 
gauge transformations as 

SA^ = d^x, (3.11) 
where x is dimensionless. The gauge transformation of is then 

^^--j^X- (3.12) 
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We have rescaled % and to fix the coefficients in these transformation rules. The 
fields and must always appear in the gauge-invariant combination 

A^^A^ + Md^cf>, (3.13) 

which can also be thought of as the gauge covariant derivative of 0: — MD^cf). 

We assume that the theory violates Lorentz invariance in the vacuum via the 
gauge-invariant order parameter 

{A^) = Mv^, (3.14) 

where is a timclikc vector of unit norm. (This defines the scale M.) In the preferred 
frame, the order parameter is (Aq). 

The coupling of the gauge field to the ghost condensate is governed by gauge 
invariance. Starting from Eq. (3.8), we can write terms 

-'-'gauged ghost 

_ FM-F^^ + £ghost(9^0 ^ AJM). (3.15) 

This does not uniquely fix the leading terms, since wc can obtain gauge invariant terms 
by replacing d^d^cp with either dfj_Au or d^A^. The general linearized Lagrangian 
expanded about the vacuum {Afj} — M5J^ is then 

-^gauged ghost = ~'~ \M'^{Aq — M)^ 

- \a,{d,Aif - la^idiAjf + \P,{dtAiY - \P2{diAof ^^'^^^ 
+ MtAodiA + HdtAof + 0{d''A^) + 0{A^). 

By the power counting of the ghost condensate effective theory, the coefficients a^, 7 
are order 1, and give rise to general kinetic terms for the vector fields. While the 
F^'^Fiiv term does not generate any independent gauge kinetic terms, for ^ 1 it 
gives the dominant kinetic contribution to the transverse modes. In the limit fl' — > 0, 
we recover the ghost condensate theory. 

Eq. (3.16) contains a "mass" term for ^0, so we can integrate out ^0 if we are 
interested in energies and momenta below the scale M. To see this explicitly, we 
compute the dispersion relation for the scalar modes. We parametrize them by 

Ao = M a, A = di(T, TT. (3.17) 

Choosing "unitary gauge" tt = 0, the quadratic Lagrangian in momentum space is 
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where 

/ (1 + g'P^)u:^k^ - g'ak' t{l - g'Ps)cok^ \ 

[ -i{l-g'(3s)u;P {I - g^[52)P + gV + g'M^ j ' ^^"'^^ 

and q; = CKi + q;2. The dispersion relation for the scalar modes is obtained by setting 
the determinant of the kinetic matrix to zero. We find two scalar modes, one with 
dispersion relation 



and one with 



J" ^-— + 0{k^). (3.21) 
7 

The first mode is gapless and corresponds roughly to the longitudinal mode a. For 

— * 

sufficiently small k we can neglect the quartic term in the dispersion relation, so the 
gapless modes travel with constant speed in the preferred frame. For g <^1, modes 
with \k\ <C gM have speed 

Cs^y/ocg^l. (3.22) 

Gapless modes with \k\ ^ gM have a quartic dispersion relation, just as the Gold- 
stone boson TT in the (ungauged) ghost condensate. In this way, this theory approaches 
the theory of Ref. [8] as (? — > 0. In the opposite limit g ^ 1 the speed is 0{1) for 
the gapless mode with <^ M, so the quartic term is never important within the 
regime of vahdity of the effective theory. 

The second mode has an energy gap of order M, and is therefore not a mode 
that is accurately described in the effective theory. We therefore write the effective 
theory below the scale M by integrating out the field Aq, which has an order- 1 overlap 
with the massive mode. We can do this without fixing any gauge, and we obtain the 
effective Lagrangian 

>Ceff = ^(M' - ^4 - i«(V • Af + 0{l/M% (3.23) 

where 

+ \^\ + 0(2. (3.24) 

9t r 9s 9^ 
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We have dropped terms in Eq. (3.23), which means that we have neglected the 

— * . 

k term in the dispersion relation for the massless scalar mode. For g <^ 1, this is 

— * . 

only justified if we restrict attention to momenta satisfying <S gM. 

What has happened is that the Goldstone mode tt has been "eaten" by the gauge 
field An and can be viewed as the longitudinal mode of the gauge field A^. However, 
unlike in the Lorentz-invariant Higgs mechanism, only gets a mass term and the 
modes parametrized by Ai remain massless. In a sense, the Lorentz- violating Higgs 
mechanism shields "electric" interactions (^o) while preserving "magnetic" ones {Ai). 
We will explore this interesting feature when we consider coupling the gauged ghost 
condensate to the standard model in Sec. 8. 



4 Coupling the Gauged Ghost Condensate to Gravity 

We now analyze the coupling of the gauged ghost condensate to gravity. We will see 
that the modification of gravity is much less dramatic than in ghost condensation. 
One aspect of this was already discussed in Sec. 3, namely that the vector modes 
mix with gravity only via derivative terms, and therefore do not modify the static 
Newtonian potential, at least in the preferred frame. Another aspect comes from 
the result of the previous section, that the gauged ghost condensate has a dispersion 
relation u"^ ~ Cgk'^ for small k, so long wavelength fluctuations are stable. Coupling to 
gravity can give only additional terms suppressed by 1/M|j, so we do not expect any 
instabilities in the presence of gravity. This is in contrast to ghost condensation, where 
the theory without gravity has a dispersion relation o;^ ~ k'^/M^, and gravity gives a 
correction Alo^ ~ — M^P/M|j that generates an instability at long wavelengths. 

Since we start with a completely generally covariant description of the theory, the 
coupling to gravity is determined simply by promoting the metric to a dynamical 
degree of freedom and replacing spacetime derivatives by covariant derivatives V^. 
For the full nonlinear theory, this is 



1 



^ (4.1) 



^The /c^ correction can be obtained by keeping the term 



which is enhanced for smah g. 
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where £eh is the Einstein-Hilbert action. 



4-1 Linear Theory 

Expanding about flat spacetime 

9^cl^ = Vnu + V' (4-2) 

the quadratic terms in the Lagrangian are 

1 

47 



>C = £eh - -j^F^^'F^. + \M\A^ - - Mf + ■ ■ ■ (4.3) 



The weak gauging of gravity affects the dispersion relation for the scalar modes only 
by terms suppressed by 1/Mpi, so we can still integrate out the massive mode as 
before. The leading quadratic terms in the effective Lagrangian for small g is then 

Aff = >Ceh + ^^i^tAif - ^F^j - + • • • ■ (4.4) 

Here 

where — a/2Mpi<I> = Mpi/ioo/v^ is the canonically normalized Newtonian poten- 
tial. Note that the (Vj^j)^ term contains mixing between Ai and hoo- 

The effective Lagrangian in Eq. (4.4) is fully invariant under gauge transformations 
and diffeomorphisms: 

SAi = da, (4-6) 

In fact, the quadratic terms in Eq. (4.4) are the most general ones invariant under 
these symmetries, and armed with the foresight to integrate out Aq, we could have 
used Eq. (4.6) as the starting point for our analysis. Indeed, this is exactly what 
we suggested in Sec. 2 when we posited the existence of a ai field. One advantage 
of thinking of Eq. (4.4) as coming from gauged ghost condensation is that the weak 
gauging procedure explains why gt and gs are nearly degenerate. 
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The modified dispersion relation for tlie scalar excitation can be easily derived. 
Parameterizing Ai = dia and choosing Newtonian gauge, we have 

where 

V 2 gMpi 

and we have taken the small g limit and dropped the subleading terms suppressed by 
g^.^ Setting the determinant of the kinetic matrix to zero, we obtain 

u'^c^k^ (4.9) 

where 

cl^a{g^-gl), g^^eg^-^^. (4.10) 

We can see that if (? > gc^ this is a healthy dispersion relation and there is no instability 
even at arbitrarily long wavelengths. If we also keep the leading contribution at 
C(A;'^/M^), then the dispersion relation is 

c^^c^P + a^, (4.11) 

which reproduces the behavior of ghost condensation for ^ gM. 

The dispersion relation of the tensor modes is also modified because (Vj^-,)^ 
contains a term ^{Ao)dthij, which contributes to the time derivative of hij. As a 
result, the speed of usual gravitational waves is modified by 0{a2M'^ /Mpi). If the 
speed of the gravitational wave is smaller than the speeds of ordinary particles, there 
is a strong constraint from cosmic rays [27, 28]. 

4-2 Modification of Gravity and Observational Constraints 

Because of the mixing between the Goldstone boson and the metric tensor, the grav- 
itational potential is modified. In the static limit, cu ^ 0, the only modification is a 
simple renormalization of Newton's constant, 

qo 2 

^^=1^ = 87rM|,(l-e2)' ^^"^2) 



^The second term in the (1,1) entry gives the leading fc^ term in the dispersion relation and 
therefore must be kept. 
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as one can see directly from Eq. (4.7). By itself, this does not lead to any measurable 
effects. 

There are non-trivial effects if the source is moving relative to the ether rest frame, 
as expected in realistic situations. For a source of mass Mq moving with velocity v 
relative to the ether rest frame, the stress-energy tensor is 

Too^p^MQS'^^\x-vt), (4.13) 

and Toi and are down by powers oi v = \v\ and can be ignored ior v <^ 1. The 
gravitational potential can be easily obtained through the propagator by inverting 
the kinetic matrix in Eq. (4.7), 

m = ill . (4.14) 

The Fourier transform of the source is 

p{(jj,k) = 2TiMQd{(jj -k-v). (4.15) 
The Newtonian potential is obtained by performing the inverse Fourier transform of 

— * 

the propagator after substituting u by k ■ v. Depending on the source velocity 
in the preferred frame v relative to the velocity of the scalar Goldstone mode = 



a{g'^ — 5(2), we get very different results. It is important to distinguish the two 
different cases 

For V < Cg, we can expand in powers of v/cs and obtain the angular-dependent 
gravitational potential 



V{r,9) = 



1 + ^ I l + —sin^ 9 + OivVct) 



(4.16) 



where cos 9 — r-i) is the angle measured from the direction of the source moving with 
respect to the ether. In the parametrized post-Newtonian (PPN) formalism [31] this 
velocity effect with respect to the preferred frame corresponds to the PPN parameter 



PPN 
"2 ' 



aPPN ^ ! = (A 17) 



The observational bound on o;?^^ is 



a 



PPN 



< 4 X 10"^ (4.18) 
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from the alignment of the solar spin axis and its ecliptic [29, 30, 31]. This provides the 
strongest bound on the scale of Lorentz symmetry breaking in this case.^ For ag"^ close 
to 1 we have M <^ 10^^ GeV, and the constraint is stronger (M <^ 10^^^/ag'^ GeV) for 
smaller g. This bound was also considered in Ref. [33] using the parametrization of 
Ref. [14] (which corresponds to a decouphng hmit of our theory as shown in App. A). 
A more systematic derivation of the modification of gravity and discussion of other 
PPN parameters can be found in Appendices B and C. In particular, we find that 
ftg is the only modification of General Relativity at post-Newtonian order, showing 
that gauged ghost condensation yielding a very mild modification of gravity. Our 
analysis also exhibits the applicable range of the PPN formalism in this theory. It is 
not enough just for v to be small because the expansion parameter is actually v/cg. It 
requires v < Cs and the mixing parameter e to be small. The above result holds for 
any theory where Lorentz invariance is broken by the vev of a vector field (such as 
Refs. [11, 12, 13, 14, 15, 16, 17, 18, 19, 20, 21, 22, 23]). 

On the other hand, ii v > Cs, the PPN expansion breaks down. We expect that 
there will be Cerenkov radiation into the Goldstone field. By causality, modifications 
of gravity happen only in a cone with an angle 6 = sm~^{cs/v) behind the source.'^ 
For V ^ Cs, the cone is narrow and we may not see modifications of gravity in the solar 
system where the most stringent bounds come from, if the cone lies outside the ecliptic 
plane. (Of course, there will be some anomalous acceleration if an astrophysical 
object happens to move into the cone of shadow of a gravitational source.) The most 
interesting effects in this case probably come from the energy loss due to the Cerenkov 
radiation. Because the coupling of the Goldstone mode to gravity is proportional to 
M, the Cerenkov radiation due to gravity scales as M^. The abihty to raise the 
Lorentz symmetry breaking scale M in gauged ghost condensation makes this effect 
interesting. In fact, it provides a significant bound on M. In contrast, in ungauged 
ghost condensation, the bound M < 100 GcV renders this effect totally irrelevant. 
However, to discuss the Cerenkov radiation from moving stars or planets, we needs 
to first understand the nonlinear effects and where they become important, so we 
postpone the discussion of Cerenkov radiation to Sec. 7 after we discuss nonlinear 
effects. 

^There is also a bound from the rate of the cosmic expansion during the Big Bang nucleosynthesis 
due to the different effective Newton's constant for the cosmological evolution [32], but it is much 

weaker. 

''This is true for distances larger than {gM)~^. For gM < \k\ < M, the fc^ term in the dispersion 
relation becomes more important than the P term, and the scalar velocity is enhanced and becomes 
k dependent, Cg ^/a\k\/M. The distance scale where this is relevant is very small if M is high and 
g is not too small. 
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5 Nonlinear Effects 



Up to now we have confined our analysis to the linear order in the fields. In this 
section, we show that nonlinear effects can be important for sufficiently strong grav- 
itational fields, and analyze the dynamics in the nonlinear regime both analytically 
and numerically. We show that the nonlinear dynamics has a simple physical inter- 
pretation as the dynamics of a charged fluid. We use this picture to give a qualitative 
understanding of the nonlinear dynamics. In particular, we argue that the caustic 
singularities found in ungauged ghost condensation are not present in the gauged 
case. 



5. 1 The Nonlinear Regime 

Recall that in ungauged ghost condensation, nonlinear effects are important for all 
interesting gravitational sources, such as stars and galaxies [10, 34]. The physical 
reason for this is that the ghost condensate gravitates like a fluid that obeys the 
equivalence principle, and therefore has a gravitational response time of order 

r 

^grav ~ -7^- (5.1) 



where r is the distance to the gravitating source. This effect appears only at nonlinear 
order; at linear order, the ghost condensate has a highly suppressed response time 
due to the a;^ = ak^ /M^ dispersion relation: 

^linear ^ — 1^ ■ (5.2) 

'a 



When fgrav ^ ^linear, the nonlinear effects dominate. These nonlinear effects were 
studied analytically and numerically in Ref. [10], and it was shown that regions with 
X — 1 < generally shrink to small size, resulting in a breakdown of the effective 
theory. This is not necessarily a disaster for the theory, since the energy involved in 
these singular regions is very small, but it is unfortunate that this behavior cannot 
be understood in the effective theory. 

In the case of gauged ghost condensation, we expect that it also gravitates on a 
timcscalc tgrav, but the linear timescale is given by the more conventional dispersion 
relation uj"^ = clk^: 



s 



T 

^linear ~ — ■ (5-3) 
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We therefore expect that gravitational effects become important when tgrav ^ ^linear, 
or 

* > (5.4) 

This means that nonhnear effects become important only for sufficiently strong gravi- 
tational fields. This will be important in subsequent sections when we discuss bounds 
on the gauged ghost condensate. 

The result Eq. (5.4) can be derived as follows. Expanding X — 1 to include the 
leading nonlinear interaction, 



the leading terms in the effective Lagrangian are 

M2 / - , _ 1 ro\^ 1 ^ „„„ a 



C = -{a,- M* - -A^) - -^F^F"^ - ^(V . (5.6) 
Integrating out the massive mode, 



we obtain 



A = M$ + — ^, (5.7) 



1 „o 1 

where we define the "electric field" 



E = dtA- MV<^-^^^^^. (5.9) 
The leading terms in the nonlinear effective Lagrangian are therefore 

^=^^^-^4-f(V-^1. (5.10) 

The leading nonlinear effects are contained in the V(^^) term in E. 

Nonlinear effects become important when the V(^^) term becomes comparable 
to the time derivative term: 

Mr 

^nonlinear (^-H) 

A 

The field amplitude A induced by a gravitational field $ is given by 

A^mV^, (5.12) 

which gives 



r 

^nonlinear ~ ~ ^infalb (5.13) 



as anticipated. 
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5.2 Fluid Picture 



We now show that the nonhnear dynamics of the effective Lagrangian Eq. (5.10) 
has a natural fluid interpretation, similarly to the one found for ungauged ghost 
condensation in Ref. [10]. This will be very useful in understanding the nonlinear 

dynamics of the theory. 

We begin by noting that in the relativistic formulation, the field defines the 
local preferred rest frame. In the effective theory with integrated out, the theory 
naturally defines a fiuid with a local 3- velocity given by 

= (5.14) 

The definition of E Eq. (5.9) can then be written as 

where 

R^d, + v-V (5.16) 

is the time derivative along the worldline of a fiuid particle (also called the convective 
or Lagrangian derivative). Eq. (5.15) has the form of Newton's law for a fiuid particle, 
with gravitational, electric, and magnetic forces on the right-hand side. 

The equations for the electric and magnetic fields arise from the Ai equation of 
motion of the effective Lagrangian, and can be written as 

dtEi + djFij = -VidjEj - ag^MdidjVj. (5.17) 

This has the form of Ampere's law with an unconventional current density on the 
right-hand side. Another difference from conventional electrodynamics is the absence 
of Gauss' law. This is expected, since the Higgs mechanism imphes that all 3 com- 
ponents of Ai represent physical degrees of freedom. 

As discussed in the previous subsection, for sufficiently strong gravity we can 
neglect the ag"^ term in Eq. (5.17), at least as long as spatial gradients are not too 
large. In this case, there is a simple solution with Ei = 0, Fij = in which the 
fiuid particles follow gravitational geodesies. This is the relevant solution in the case 
where the fiuid is initially "at rest" in the presence of a gravitating source. Just as 
in the ungauged case, the fiuid particles will "fall in" toward the source, giving a 
very direct and physical picture of why the gravitational time scale is relevant for the 
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nonlinear dynamics. The subsequent evolution of the fluid will in general give rise 
to caustic singularities. Near these caustic singularities, the higher-derivative ag^ is 
important, and may resolve the singularity. In the fluid picture, the "electric" and 
"magnetic" forces are becoming important near the would-be caustic, and may cause 
the incoming fluid particles to "bounce." This is the question we address next. 



5.3 Would-be Caustics 

To address the question of caustic singularities in the nonlinear dynamics, we restrict 
attention to situations with spherical, cylindrical, or planar symmetry. We can treat 
these simultaneously by using the spatial metric 

^ ^^2 ^ ^2^^2 ^ J2 dxl, (5.18) 

p=l 

where 

{0 s = (planar symmetry) 

dO"^ s = 1 (cylindrical symmetry) (5.19) 

d6^ + sin^ Odip"^ s — 2 (spherical symmetry) 

Since only Ar is nonzero, we have Fij — 0, and "Newton's law" simphfles to 

^ = -dr^ - ^Er, (5.20) 

Dt M ' ^ ^ 

where D/Dt — do+Vrdr- "Ampere's law" Eq. (5.17) can then be written in the form 

^^{r'Er) = -r^jr, (5.21) 



where 



ag^Mdr 



—dr{r^Vr) 



(5.22) 



Eq. (5.21) says that the comoving "flux" r'^E^ is changing according to the "current 
density" jr- In particular, when is negligible, the "flux" is conserved, a fact that is 
useful in interpreting the numerical results below. 

We now consider an initial condition that can lead to a caustic singularity. That 
is, we assume that initially Vr < near r — 0, corresponding to a situation where the 
fluid particles are heading toward the origin. More precisely, we assume that near 
r = 0, we have 

-cir + C2r^ + 0{r^), (5.23) 
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with ci > 0. Because we expect the hnear growth of Vr away from r = to degrease 
with r, we also assume C2 > 0. In this case, we find that near r — 

jr = (s + 2)ag^Mc2 > 0. (5.24) 

Therefore, Eq. (5.21) imphes that the "electric field" Ej. tends to decrease with time. 
In particular, if E'^ = initially {e.g. if the fluid is "at rest") then later Ej. < 0. 

The sign of Ej. is crucial for the fate of the would-be caustic. For Ej. < 0, the 
"electric" force in Eq. (5.20) is repulsive and can cause the incoming particles to 
bounce. To see the condition for this, we rewrite Eq. (5.20) as a function of the 
comoving radial coordinate r as 



- = -- 



= -TjEr. (5.25) 



where we have neglected the gravitational "force." We therefore have 

This can be interpreted as the work-energy theorem for fluid particles. If E'r < and 
the "work" integral on the right-hand side diverges as r/ ^ 0, then the particles will 
bounce. This occurs for s = 1,2 (cylindrical or spherical symmetry), since E^ ~ r~* 
neglecting the "current" contribution in Eq. (5.21). Physically, it requires an infinite 
amount of energy to compress a charged sphere or cylinder to zero size. For s = 
(planar symmetry), the "work" integral is finite, corresponding to the fact that a 
finite amount of energy is sufficient to compress a plane of charge to zero size. Indeed 
numerical simulations (discussed below) confirm that caustic singularities do form for 
planar symmetry, but not for cylindrical or spherical symmetry. However, we expect 
that departures from perfect planar symmetry will be important near the would-be 
caustic. Specifically, small fluctuations of the plane symmetric collapse will grow, 
and we expect the layer to fragment before a planar caustics occurs. (See App. D 
for a perturbative analysis supporting this picture.) After the fragmentation, lower- 
dimensional caustics do not occur by the argument above. 

The fact that the "current density" j,- generates a negative E^ depends crucially 
on the non- vanishing second derivative of at r = 0. In fact, the nonlinear equations 
have exact solutions 

A - E - ^J}^Z^ (527) 

where = 1 or 2/(s + 1) and is an arbitrary constant. For these solutions jV = 0, 
and nothing prevents the solutions from collapsing at t = tc- These solutions arc 
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of course unrealistic because of the boundary conditions at r — >^ oo, but we might 
worry that they are approximately valid near the origin. However, we expect that 
the scaling solutions Eq. (5.27) are dynamically avoided. This is because the second 
derivative of Ar must be negative away from r — 0, and this propagates to the origin, 
making > and preventing collapse. 

5.4 Numerical Results 

We now turn to numerical investigations of the nonlinear equations. We begin with 
the plane symmetric case (s = 0) with no external gravitational potential = 0). 
For simplicity we assume the symmetry under reflection r —r. In this case A^ and 
Er are odd functions of r. As the initial condition at t — we set 

Ar\t=o = 2nre-'-\ Er\t=o = 0, (5.28) 

where ^2 is a positive constant. By performing numerical simulations, we find that for 
some values of ag^ and fl the would-be caustic bounces, but much later it recoUapses 
and forms a caustic that exits the realm of vahdity of the effective field theory before 
it bounces. This is shown in Fig. 1. For other values of ag"^ and Q, even the first 
caustics at i ~ does not bounce until the system exits the regime of validity of 
the effective field theory. This confirms the conclusion of the analytical argument 
above. As discussed there, we do not expect planar caustics to form in the realistic 
case. 

We now turn to systems with cylindrical and spherical symmetry. Figs. 2 and 3 
show results of numerical simulation in cylindrical and spherical cases with $ = 0, 
respectively. It is evident that the excitation near the origin diffuses and the system 
asymptotically approaches the trivial solution Ar,Er = 0. These results confirm our 
expectation that caustics with cylindrical or spherical symmetry do not occur. 

So far, we have been working in the absence of gravity. Let us now consider effects 
of an external gravitational potential. We still neglect gravity generated by excitations 
of the gauged ghost condensate. For simplicity we assume spherical symmetry (s = 2). 
A typical numerical result is shown in Figs. 4 and 5. We see that in the final stationary 
configuration, the repulsive "electric force" E^ cancels the attractive gravitational 
force {i.e. A^ 0, E^. ^ -dr^). 

6 Black hole accretion 

In this section we analyze perturbations of Schwarzschild geometry and estimate the 
rate of mass increase of the black hole due to accretion of the gauged ghost condensate. 
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Fig. 1. A result of numerical simulation for planar symmetry (s = 0). The variable 
Ar is plotted as a function of r and t. After an initial bounce, a caustic singularity 
forms at late times. Note that in a realistic situation, we do not expect a system to 
exhibit planar symmetry. 




Fig. 2. A result of numerical simulation in cylindrial symmetry (s = 1). The variable 
Ar is plotted as a function of r and t. After an initial bounce, no caustics form at 
late times. 
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Fig. 3. A result of numerical simulation in spherical symmetry (s = 2). The variable 
Ar is plotted as a function of r and t. Like the case of cylindrical symmetry, there 
are no late time caustics. 




Fig. 4. A result of numerical simulation in spherical symmetry (s = 2) with an 
external gravitational potential. The variable is plotted as a function of r and t. 
Again, we observe the absence of late time caustics. 
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Fig. 5. A result of numerical simulation in spherical symmetry (s = 2) with an 
external gravitational potential. The variable Er is plotted as a function of r and t. 
At late times, the "electric field" cancels the gravitational force, halting the accretion 
of the gauged ghost condensate and protecting the system from caustics. 
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extending the analysis of Ref. [35] for the ungauged ghost condensate. 

Before starting the discussion of black holes, let us first consider a spherically 
symmetric, static star surrounded by a gauged ghost condensate. In the previous sec- 
tion we have seen that the gauged ghost condensate with spherical symmetry quickly 
approaches a regular, stationary configuration with or without an external gravita- 
tional potential. In the fluid picture, each fiducial fluid particle follows "Newton's 
law" and approaches the center within the dynamical time set by the external gravi- 
tational potential. Near the center, the "electric field" Er builds up and the infalling 
flow of fluid particles bounces. Part of outgoing flow produced by the bounce may 
come back to the center at late time but bounces again. After several (if not one) 
bounces, the system settles down to a regular, stationary configuration. Inclusion of 
gravitational backreaction just induces the "renormalization" of Newton's constant 
(see Eq. (4.12) for the renormalization in the linear theory) and does not change this 
qualitative behavior towards the stationary configuration since the stress-energy ten- 
sor of the final stationary configuration is time-independent and does not have any 
time-space components. Note that for the ghost fluid around a star to settle down 
to the stationary configuration, the existence of a regular center is crucial: bounces 
happen because of the boundary condition at the center. 

The ghost condensate surrounding a black hole behaves very differently from that 
around a star. This is because the boundary condition at the horizon of a black hole 
is completely different from that at the center of a star. Fiducial fluid particles just 
go through the horizon and continue to infall. The flow of ghost fluid should still 
be regular from the viewpoint of infalling comoving observers: shell-crossing type 
singularities are avoided by bounces due to temporal build-up of the local "electric 
field" . Thus, if we neglect the gravitational backreaction, the ghost condensate sur- 
rounding a black hole should approach a congruence of free-falling geodesies in the 
fluid picture and this should happen roughly within the dynamical (or Kepler) time. 
The stress energy tensor vanishes if a = and is suppressed by the factor aM^/Mp^. 
Hence, the gravitational backreaction does not change this qualitative behavior hap- 
pening within the dynamical time. Thus, the congruence of free-falling geodesies is 
a very good approximation to the full solution including gravitational backreaction. 
However, if we are interested in physics of much longer time scales, the gravitational 
backreaction may build up to appreciable amount. Therefore, we would like to esti- 
mate the mass increase of a black hole due to accretion of the ghost condensate. In 
the following, we first find the approximate solution representing the congruence of 
free-falling geodesies and then estimate the gravitational backreaction by perturbative 
expansion with respect to aM^/Mp^. 
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In order to investigate the behavior of the gauged ghost condensate in a black hole 
background, we need to choose the time variable carefully. For example, if we use the 
background Killing time as a time variable then there would appear infinite blue-shift 
on a black hole horizon. With the infinite blue-shift, the local physical energy scale 
near the horizon would easily go beyond the scale M. In particular, we would not be 
allowed to integrate out the field Aq, which has an order- 1 overlap with the massive 
mode in Eq. (3.21). For this reason we should use a different time coordinate for 
which there is no large blue-shift. In the following analysis we shall use an infalling 
Gaussian normal coordinate system, in which there is no blue-shift since the time- 
time component of the metric is —1 everywhere by definition. In this coordinate 
we can safely use the effective theory with Aq integrated out. As explained in the 
previous paragraph, the gauged ghost condensate surrounding a black hole quickly 
settles down to a congruence of infalling geodesies in the fiuid picture. In other words, 
the solution in the zeroth order in aM'^/Mpi is given hy (j) — M'^r, where r is the 
time coordinate in the infalling Gaussian normal coordinate system. 

The unperturbed Schwarzschild metric in the infalling Gaussian normal coordinate 
called the Lemaitre reference frame [36] is 

ds'' = -dr^ + ^^^dx' + r^T, x)d^l\ (6.1) 
r(T, x) 



where 



r(T, x) — 2mo 



3 / r 

X 



2/3 



(6.2) 



.2 V 2mQ, 

There is nothing bad on the future (black hole) horizon and the coordinate system 
covers everywhere in the region v > — oo, where v is the ingoing Eddington-Finkelstein 
null coordinate [37] . The metric becomes ill only on the curvature (physical) singular- 
ity at r = 2mQX. We consider spherically-symmetric, time-dependent perturbations 
of the Schwarzschild geometry. We shall still use an infalling Gaussian normal coor- 
dinate system: 

ds"" = -dr"" + ^ ^ , dx-" + e2^^(-'-)r2(r, x)dn\ (6.3) 

r(T, x) 

where r(r, x) is given by (6.2). We consider 5i and 82 as perturbations. Up to the 
linearized level, the Misner-Sharp energy [38] is expanded as 



Mms - " " 



d^{e^^r)d^{e^^r) = M'^i{mo + rm) , (6.4) 



where 
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Note that the Misner- Sharp energy agrees with the ADM energy at large distances 
and, thus, measures the strength of gravity. In the following we shall eliminate 5i 
from all equations by using this expression and consider mi and S2 as dynamical 
variables. As for the gauged ghost condensate, we introduce linear perturbation 
A{t,x) as Aidx^ — A{T,x)dx. Einstein's equation gives a set of coupled partial 
differential equations for the three variables mi, ^2 and A. Fortunately, it is possible to 
decouple S2 from equations for mi and A by taking the following linear combinations 
of the linearized Einstein tensor: 



{2mof 



2mo 



2mo 

1/2 



3/2 r 



G 



a 



Lx 



±± 



(2mo)' 



2mo 



mo 

drTTli 

mo 



(6.6) 



By applying the formula in Appendix B to 

N=l, /3 = 0, 

(Om '\3„25i(r,a;) 

qi.dx'dx^ = ^ "^Z ^ , dx^ + e2^^(^'^V2(T, x)dVL\ 

r{T,x) 

Aidx^ = A{t, x)dx, 



(6.7) 



we can obtain the stress energy tensor of the gauged ghost condensate. 

In this setup, we expect that accretion of condensate to a black hole is due to 
the a term. Hence, we perform a perturbation analysis with respect to aM"^ /Mpi, 
considering mi, ^2 and A as quantities of the first order in the perturbation. At zeroth 
order in aM'^/Mpi, all components of the stress energy tensor vanish and the Einstein 
equation is automatically satisfied. This means that there is no mass increase when 
aM'^/Mli = 0. 

At first order in aM'^ /Mpi, the non- vanishing components of the stress energy 
tensor are 



T 



±± 



T 



±x 



M2 


1 r 


(2mo)2 


g"^ 2mo 




8mo V 






CD 


8(2mo)4 




9aM^ 





-dAA+ — 
2g^ 



2mn\^/2 



drA+ — 



9a /2moy 



8(2mo)4 \ r J 



(6.8) 
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Hence, we obtain 



T 



2mo (2mo)' 



1 ''.a,a..4 + Af^V"a.^-''"^"''" 



(yf^ 2mo ^ 2^1^ \ r J ^ 8 V 



r 



+ (2mo)5 8(2mo)n r J I mo;' ^ ^ 

Thus, Einstein equation implies that 

1 3^^+ 5 /5«5y'^g^^^9a/5^^'_A^/2^^ ^^^^^^ 



g"^ 2mo 2gf2 V r / 8 V r y V r / mo 

and 

The first equation should be considered as an equation for A for a given mi. The 
second equation does not include A and, thus, can be easily solved w.r.t. mi. The 
general solution for mi is 



mi 9aM2 



mo 4M' 



2 

Pi 



1 , / r 
+ - In 



2mn 2 V 2mi 



(6.12) 



where C'(a;) is an arbitrary function. Note that r parameterizes each geodesic as 
the proper time measured by a comoving observer and x parameterizes the congru- 
ence of geodesies. Thus, the function C{x) represents difference between evolution 
along different infalling geodesies and thus the initial condition on an initial spacelike 
hypersurface. 

This result is equivalent to the corresponding expression (22) in [35] for the un- 
gauged ghost condensate. To see the equivalence, notice that x here corresponds to 
x+ in [35]. Thus, the accretion rate up to first order in aM'^/Mpi is the same as that 
in the ungauged case. In particular, as shown in [35], the leading late-time behavior 
is 

mo 4M|,, V4mo/ ^ ^ 

where v is the advanced null time coordinate normalized at infinity (the ingoing 
Eddington-Finkelstein null coordinate). To obtain this formula, the asymptotic be- 
havior of the function C{x) for large x has been determined by the assumption that 
the initial value of mi on an initial spacelike hypersurface does not diverge at large r, 
and the limit v ^ r has been taken. This formula is valid in the regime where mi is 
sufficiently small compared to mo. In order to see the behavior beyond this regime. 
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let us "renormalize" the parameter mo- For this purpose we note that 



1 djM'jH) ^f3\'/'aM^ 
^4/3 ^(^2/3) -'UJ Ml, ^^-'^^ 

for the black hole mass Mbh = (^o + mi)Mpi. This gives the solution 

mTh ~ Ml'' + 2 (5) "-^(Ml,vf^ (6.15) 

for the mass Mbh surrounded by the sphere at the radius r at the time where Mq 
is the initial value of Mbh at = 0. This formula is vahd if the condition i> » r is 
satisfied. We shall check this condition when we apply this formula to a stellar-mass 
black hole below. 

Convincing evidences for stellar-mass black holes are provided by X-ray binaries, 
such as Cygnus XI. Some candidates for the stellar-mass black holes are listed in 
Table 1 of [39]. On the other hand, ages of those stellar- mass BHs are less certain. 
For example, XTE Jl 118+480 is thought to be a black hole with mass Mbh ~ 
and binary separation r ~ 3i?0, but its age is estimated to be either ~ 240 Myr 
or ~ 7 Gyr [40], depending on whether it was kicked by a supernova explosion or 
was ejected from a globular cluster. (GRO J1655-40 has Mbh ~ 5Mq and its age is 
estimated to be ~ 0.7 Myr [41]. This would give a weaker bound.) 

The observation of XTE Jl 118-1-480 imphes that the mass measured at r ~ ?>Rq 
is Mbh ^ 7Mq at t ~ 240 Myr or 7 Gyr, while the theory of stellar evolution 
says that Mq must be larger than ~ 3Mq. Note that the estimate of t should not 
change significantly even if Mbh actually evolved from ~ 3Mq to ~ 7Mq during its 
journey. We can apply the formula (6.15) to this system since M^iT /Mbh ~ 10^ and 
Milt /Mbh ^ 10^° and thus the condition V ~ t 3> r is satisfied. Therefore, we obtain 

(7Me)^/^ - {SMef > ^{Mli x 240 Myrf/\ (6.16) 

This gives the bound ^ on M as 

y/aM < 10^2 GeV. (6.17) 



^Because of the motion of the binary system, the relative rotation of the binary components and 
the non-zero spin of the black hole, the spherical symmetry assumed in the derivation of the formula 
(6.15) is violated. It is expected that these factors will increase the accretion rate. Thus, while 
(6.17) is valid as an upper bound on sJaM, it is worthwhile seeking more stringent bound by taking 
them into account. 
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7 Cerenkov Radiation from Gravitating Sources 



For a gravitational source moving faster than the sound velocity of the ether {v > Cg) 
ether Cerenkov radiation will be emitted. As we will see, the energy loss per unit 
time is at most of order A'Pv^, where M is the scale of the gauge ghost condensate, 
independently of the size of the source. In the ungauged ghost condensate, the same 
result holds [9], but in the gauged case M can be much larger, which give rise to 
interesting phenomena and constraints from this effect. 

Consider a non-relativistic source 



V2Mpi 

As shown in App. E, the energy loss rate can be calculated as 

1 



\/2M] 



-p{x,t)^%x,t). 



PI 



At the linearized level, 

^"(u, k) 



1 aM^ 



(7.1) 



(7.2) 



(7.3) 



where we have ignored the k^ piece in the denominator of the second term in the 
bracket because it is never important for the large astrophysical bodies we consider 
here. 

For a star of mass Mq moving on a straight line with velocity v: 

p{uj, k) = 2t:Mq5{uj - k ■ v)f{k), (7.4) 

where f{k) is the dimensionless form factor given by the Fourier transform of the 
mass density distribution. For example, a uniform-density sphere of radius Rq gives 

3fsin /ci?0 — kRQ cos /ci?©) 



Substituting Eq. (7.4) into Eqs. (7.2) and (7.3), we have 
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(7.5) 



\m\', (7.6) 



'PI (k ■ V + isY — c^k^ _ 

where we have used the ie prescription for the retarded Green's function. Note that 
the integral is only nonzero due to the poles of the second term, i.e. when the prop- 
agator is on shell. This occurs only for v > Cs, and we obtain 



E = 

-•^src 



IGttM^ Jo 



kdk\f{kf 



(7.7) 
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for u 3> C5. For a source of size Rq, we expect 



f 

Jo 



kdk\f{kf\r.-^. (7.8) 



However, as we discussed in the previous section, nonlinear effects can be impor- 
tant near a gravitational source, and may cut off the integral Eq. (7.7) at radius larger 
than Rq. Using the estimates of the previous section, we see that nonlinear effects 
become important within a distance 

r <i -Rnonlinear ~ 2 ' C^"^) 

In the ungauged ghost condensate, the parametric formula for the rate of energy loss 
in the linear regime continues to be valid in the nonlinear regime [10], and we find 
the same result here for essentially the same reason. Consider a fluid particle moving 
past a source with speed v and impact parameter r. In the impulse approximation, 
a fluid particle gets a kick perpendicular to its initial velocity due to gravity 

Here we have neglected the "electric" force on the particle. This can be estimated 
using 

J-E^ ^ ^^At ^ ^ ^W. (7.11) 

We see that neglecting the "electric force" is justified for v > Cg. The impulse approx- 
imation is valid as long as the perpendicular distance travelled by the fluid particle 
is less than r, which gives 

. & fl*.. ~ (7.12) 

Note that it!drag < -^nonlinear iOT V > Cg. FOT T < i^dragj We expect the fluid tO 

be "dragged" with the source, and no Cerenkov radiation will be emitted from the 
dragged region. 

The energy loss induces an anomalous acceleration in the direction of the ether 
wind, 

Esrc aKMQM"^ 



Because it is proportional to Mq/Rq, the Sun gets the largest anomalous accelera- 
tion in the solar system. These accelerations will modify the orbits of the planets. 
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In particular, they will cause misalignments of the orbital planes if the anomalous 
accelerations are out of the ecliptic plane. For comparison, the typical relative accel- 
erations between the planets and the Sun are 10"^^ — 10~^^ GeV, the Viking ranging 
data constrains any anomalous radial acceleration acting on Earth or Mars to be 
smaller than ~ 10~^^ GeV while the Pioneer anomaly corresponds to apioneer ~ 10~^^ 
GeV [42] . If we require the anomalous accelerations to be smaller than ~ 10""^^ GeV, 
we obtain an upper bound on M: 

M ^ 10^°(afi;)-i/2 GeV. (7.14) 

Another potential constraint outside the solar system comes from the period 
change of binary pulsar systems. Note that even if the velocities of the pulsars are 
smaller than Cs, the Goldstonc boson can still be emitted through the quadrupole 
radiation. To simplify the problem, we assume that two pulsars of equal mass Mq 
move in a circular orbit separate by a distance 2ro with an angular frequency Uq. The 
angular frequency and the orbital velocity vq are related by vq — ouoTo- We also ignore 
the velocity of the center of the system. 

The energy loss due to the multipole radiation in the case where Vq < Cg or 
Cerenkov radiation in the case where vq > Cg can be calculated in the same way 
described in App. E. Analytic formulae can be obtained by taking different limits of 
the Bessel functions for vq <^ Cg or ^ c^. 

For Vq — ujqTq <tC Csi the energy loss is dominated by the quadrupole radiation. As 
derived in App. E, it is given by Eq. (E.13). 

15nM^,ricl~ 15 cl ' ^ ^ ^ 

On the other hand, the energy loss due to usual gravitational waves is given by [31] 
_ ^ hnn\- 2Gi.Mo5_ 2^^7rMg,^o^° 



'PI 



r, 



where Qi/s are the quadrupole moments of the gravitational source. Taking the ratio, 
we have 



PPN 




(7.17) 



The observed period change of the binary pulsars PSR 1913+16 is consistent with the 
energy loss entirely due to the gravitational wave radiation [31], so the above ratio 
must be smaller than 1. Given their orbital velocity vq ~ 10~^, we can see that for 
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vo -C Cs (say, > 10 fo), this does not give a stronger bound than the a^^^ bound 
from the solar system. For vq getting close to this can compete with the solar 
system bound, but our approximation starts to break down. 

Next, let us consider the opposite limit, Cg <^ vq. The energy loss is given approx- 
imately by (see Eq. (E.15) of App. E) 



■ aM'M^ ^ 2nvo 



J, / 2nvQ \ 



(7.18) 



The problem is that the arguments in the form factor take values much bigger than 
1/ro in this case. As we discussed before, nonlinear effects are already important 
at that length scale and the region of the ether being dragged by the pulsars is 
approximately of size Tq. The effective form factor must be very suppressed and is 
difficult to calculate from first principles. Physically, the binary system emits ether 
waves with angular frequencies that are integer multiples of 2u!o- Because the sound 
velocity of the ether is much smaller than the orbital velocity, the wavelengths of the 
ether waves are much shorter than the size of the system. Therefore, the amplitudes of 
the ether waves must be highly suppressed. If we assume a power law suppression that 
f{2nvQ/ CsTq) ~ (2nfo/cs)~P with p > 1, then the sum in Eq. (7.18) is ~ (cs/vq)^^"^, 
and the energy loss would be 

F "^^'^^'^ r V'"' 1 ( \ (^^''\ (^^\ 



AT:M^^rlcs\vQj VlO^o GeVy V vq ) \vqJ dt 

(7.19) 

Obviously for Cg <^ Vq this does not give a significant constraint. For Cg getting 
close to Vq it may become a competitive constraint depending on the power of the 
suppression factor, but again our approximation starts break down. 



8 Couplings to the Standard Model 

Given that gauged ghost condensation yields such a modest modification of gravity, 
the strongest bounds on this model of Lorentz breaking could come from direct cou- 
plings between Ajj, and the standard model. As we will see, most couplings could 
be forbidden by discrete symmetries, but we will consider what would happen if the 
gauged ghost condensate were not confined to a hidden sector. 

Even in the absence of gravitational couplings, the Lagrangian in Eq. (3.23) is 
fascinating because it is the effective field theory description of an Abelian gauge 
theory in a Lorentz-violating Higgs phase. To our knowledge, the dynamics of this 
phase has not been studied in the literature. Like the Lorentz-invariant Higgs phase. 
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the gauge boson eats a Goldstone and therefore has three polarizations, but what is 
especially bizarre about this system is that all three polarizations are massless. There 
are long-range "magnetic" interactions in this theory but "electric" interactions are 
exponentially suppressed. In other words, in order to produce a healthy Newton's law 
for gravity, we had to sacrifice Coulomb's law for this ghost-electromagnetic U{1). 

Note that this U (l)ghost theory is different from other Lorentz- violating U{1) gauge 
theories considered in the literature. In theories based on Ref. [1], the Lagrangian 
contains gauge invariant but Lorentz- violating terms like kaP'ysF'^'^F^^ , where kai3^s 
is some fixed Lorentz-violating tensor. Because gauge symmetry is maintained, the 
equations of motion for electric and magnetic fields are modified without introducing 
new propagating degrees of freedom. A model that does violate gauge symmetry 
is given in Ref. [43] which considers a "magnetic" Higgs phase where Coulomb's 
law is virtually unchanged. Starting with the proposal of Ref. [44], there has also 
been speculation that the electromagnetic field could arise as Goldstone bosons from 
spontaneous Lorentz breaking. In contrast to these other theories, gauged ghost 
condensation is a model where a U{1) gauge field enters an "electric" Higgs phase 
triggered by a Lorentz-violating vev for a charged scalar field. 



8.1 Catalog of Allowed Couplings 

At minimum, even if we do not include direct interactions between standard model 
fields and Afj,, we expect graviton loops to generate interactions of the form 

AC ^ c^T^^A.A., (8.1) 

where T^,^ is a symmetric dimension four standard model operator, and c is expected 
to be 0{1). Note that there is no constraint from gauge invariant on this couphng 
because A/j. — + Md^cf) is a C/(l)ghost gauge invariant combination. The l/Mpj 
suppression arises because A^ contains d^4> so this interaction is actually dimension 
eight, and the factors of M are inserted because A^j, is defined with mass dimension 
Setting A^ to its vacuum expectation value (^o) = 

AC ~ 0^7*^°. (8.2) 

If we take T'^" to be the stress-energy tensor for a standard model fermion, the 
inclusion of T*^° allows the maximal attainable velocity of the fermion to differ from 
the speed of light: 

5v, ^ ^. (8.3) 
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Constraints from high precision spectroscopy and the absence of vacuum Cerenkov 
radiation bound Sv^ < 10~^^ — 10^^^ [45]. This hmits the scale of spontaneous 
Lorentz violation to be M ^ 10^^ GeV. Remarkably, for ~ 1/10 this bound on M 
is comparable to the bound from the gravitational sector in Sec. 4. 

If there is a — > —An symmetry, then Eq. (8.1) is the leading couphng between 
the gauged ghost condensate and the standard model. Relaxing this symmetry allows 
couplings of the form 

AC = qJ^A^,, (8.4) 

where J'^ is a standard model current and q is the coupling constant. In fact, this is 
just equivalent to charging standard model fields under the U (l)ghost ga-uge symmetry. 
To see this, consider the interaction 

?V AV'T'^V' = 9v i^t^ + Mdn(f))'4;YiP, (8.5) 

where i/j is a standard model fermion. If the current i/j^'^i/j is conserved, its interaction 
with dn4> can be removed by a field redefinition 

iP^iP' ^ e-^^^^'t'iP (8.6) 

into the kinetic term. In terms of the redefined field the interaction Eq. (8.5) 
becomes 

q^A^^Y^' (8.7) 

which is just a t/(l)ghost gauge interaction of a fermion with charge q^. The point is 
that exp{—iM(f)) carries a unit charge of the C/(l)ghost gauge symmetry, so the charge 
of a field under [/(l)ghost can be altered by multiplying it by powers of exp(— iM0). 
It is most convenient, however, to work in the ip' basis where the 9^0 coupling is 
removed. 

New long-range forces other than electromagnetism and gravity are strongly con- 
strained. In gauged ghost condensation the Aq field is massive, so there is no 1/r 
potential between C/(l)ghost charges in the preferred rest frame. Nonetheless, Ai is 
still massless, giving rise to additional magnetic forces. Effectively, this gives a tree- 
level contribution to the magnetic moment of a particle charged under the f/(l)ghost- 
If we assume that the ?7(l)ghost charges of the standard model fields are proportional 
to their C/(1)em charges, the correction to the Lande g factor is 



~2 

Qe '-^ghost 





q;em 



(8.8) 
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where ctghost = 5'^/(47r). Note that in theories with a normal para-t/(l) gauge group 
that "regauges" electromagnetism [46], this correction is absent, because the para- 
U{1) also shifts the electrostatic charges of the fermions. The most accurately mea- 
sured gi — 2 is for the electron with an uncertainty at the 10~^ level [47, 48]. However, 
since it serves as a definition of ckem, to constrain q^g we need the next most precise 
determination of ckem- A measurement of a^u using the atom interferometry with 
the laser-cooled cesium atoms has reached the 10"^ uncertainty level and the result 
is in agreement with the g — 2 measurement [49] . This implies a constraint 

S^^^lO-s ^ Qeg ^10-'. (8.9) 
ckem 

As we will see, this bound is much weaker than the bound from the long-range force 
derived in the next subsection: for sources moving with respect to the preferred 
frame, there is a direction-dependent pseudo- Coulomb potential. If the C/(l)ghost 
charges are not proportional to the electromagnetic charges, then the measurement 
of {g — 2) would depend on which species of fermion was responsible for setting up the 
background magnetic field, though in practice, almost all laboratory magnetic fields 
are in some way established by electrons. 

One can also consider the couplings of standard model antisymmetric tensors to 
the field tensor of the f/(l)ghost- At the lowest order, we can write down a kinetic 
mixing between the U{1) ghost and U{1)em, 

AjC = CoF^^FZ,. (8.10) 

Such a term can be removed by a redefinition of the U{1) gauge fields, resulting 
in a shift of the ^7(l)ghost charges of the standard model fields, so the effect is the 
same as discussed in the previous paragraph. Alternatively, we could imagine dipole 
interactions 

A£=^^c*fr^^^*F^.. (8.11) 

These couplings were studied in a Lorentz invariant setting in Rcf. [50], where the 
bound on A for the electron was expressed in terms of the Higgs vev {h) = 174 GeV: 

A>P-W,10^-10^GeV. (8.12) 

Ce {h) 

where the more stringent bound comes if Eq. (8.11) introduces flavor mixing between 
electrons and muons. It would be interesting to study these couplings more carefully 
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in the Lorentz- violating Higgs phase, but we do not expect significant departures from 
these bounds. 

Because we will find strong bounds on from direction dependent pseudo- 
Coulomb potentials, we want to forbid the coupling in Eq. (8.5) while still allow- 
ing new couplings other than just Eq. (8.1). One way to forbid a coupling to the 
electromagnetic current is to choose to be parity-odd, or equivalently 

P : Ao^ -A, Ai ^ Ai- (8.13) 

In the standard model, electroweak interactions violate parity, and because we expect 
the scale M to be much higher than the electroweak scale, at best wc can say that the 
coupling of A/j, to parity-even currents should be suppressed relative to parity-odd 
currents. For a fermion ^ with Dirac mass nio, we can construct the vector and axial 
currents: 

= ^r^^*, 4 = ^7m7'*- (8.14) 

By looking at the relevant triangle diagram, we estimate that the coupling A'^J/^ 
should be suppressed relative to the coupling A^J^ by m|,G'F/167r^, where G^? is 
Fermi's constant. For an electron, this suppression is 10~^^, so to a good approxima- 
tion we can ignore couplings to parity-even currents. 

With the above caveats, the leading coupling of Afj, to any light Dirac fermion, in 
particular electrons and nucleons, is 

Ant = Y^'Jl (8.15) 

where we parametrize the coefficient as a ratio of two scales M and F to compare 
with the result of ungauged ghost condensation. Setting A^ to its vev: 

This Lorentz- and CPT- violating term is the same as in ghost condensation [8, 9], 
and gives rise to different dispersion relations for the left- and right-handed fermions 
[51, 52]. In a frame boosted with respect to the preferred frame, Eq. (8.16) contains 
the interaction 

/iS-^4arth, (8.17) 

where we have identified '3>77^^I' with non-rclativistic spin density s. Assuming the 
preferred rest frame is aligned with the CMBR, we take ItTcarthl ~ 10^^, and experi- 
mental bounds on velocity-dependent spin couplings to electrons is ~ 10~^^ GeV 
[53] and to nucleons ~ lO'^^ GeV [54, 55]. 
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After integrating out ^o, we are left with the interaction 

M 

Ant^^s-.A (8.18) 

r 

The (7—^0 hmit of this couphng was discussed in Ref. [9] and it led to two Lorcntz- 
violating dynamical effects: "ether" Cerenkov radiation and a long-range spin-dependent 
potential. These effects in this model will be discussed in the next subsection. 



8.2 Dynamics of the Lorentz- Violating Higgs Phase 

The Lagrangian in Eq. (3.16) is the most general Lagrangian with an An —An 
symmetry that is spontaneously broken by the vacuum (^0) — M. If we only enforce 
Eq. (8.13), then we can add additional couplings to the gauged ghost Lagrangian such 
as 

AC ~ M{Ao - M)diA\ (8.19) 

When we integrate out Aq, however, this term docs not change the basic form of 
Eq. (3.23), and the effective Lagrangian for Ai at distances large compared to 1/M 
is 

Aff = ^2i9tAr - ^4 - i«(V • Ar, (8.20) 

where we have absorbed the effect of Pi and a2 into the parameters g and Cy, and 
has the interpretation of the velocity of the transverse vector modes. 

Going to Ai = Ai gauge, it is straightforward to calculate the AiAj propagator, 

where Cg = Cyg\foL is the velocity of the longitudinal mode of Ai. We will always 
assume that ^ c.y. The first thing to check is what kind of pseudo-Coulomb 
potential is generated between two particles carrying ^7(l)ghost charges in a frame 
moving with respect to the preferred frame. Consider a source of charge Q under 
C/(l)ghost moving with velocity v < Cs- 

J° = (9(5(3) _ jO ^ 2'KQb{u -k-v) 

(8.22) 

J ^ Qv5^^\x - vt) J = 27iQv6{u-k-v) 
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where there are corrections at 0{v'^). To this order, wc can use the a; — >■ hmit of 
Eq. (8.21). Taking the Fourier transform of the propagator, the potential between Q 
and a test charge q a comoving distance r away is 

V{r, e) = -3^^ ( 1 + cos^ e + ^ sin^ o] , (8.23) 
Svrr y g^a J 

where cos^ — f ■ v. In the hmit <^ c^,, the third term dominates, yielding an 
effective anomalous gauge coupling 

5anom ^ (8.24) 

and angular dependence with respect to the velocity of the preferred frame. Bounds 
on new gauge interactions are usually derived by noting the absence of any long- 
range forces other than electromagnetism and gravity. In particular, precision tests 
of gravity using materials with different baryon-to-lepton ratios place limits on a 
possible U{1)b-l gauge coupling of Qb-l < 10~^^ [56]. If searches for U{1)b-l gauge 
couplings are sensitive to such velocity effects, then, qganom ^ 10~^^. Assuming the 
earth is moving with velocity \v\ ~ 10~^ with respect to the preferred rest frame, then 

< lO-^'^V^, (8.25) 

where ip represents an ordinary particle p, n, or e. For reasonable values of a and 
g, this bound is stronger than Eq. (8.9), indicating that couplings between and 
vector currents must be strongly suppressed. Still, it would be interesting to search 
for the angular dependent Coulomb potential in Eq. (8.23) to see whether the bound 
on q^ could be improved. 

Even if we forbid the coupling to vector currents using Eq. (8.13), axial couplings 
are still allowed. The couphng in Eq. (8.15) between A and spin is familiar from 
ghost condensation and gives rise to a spin-dependent inverse-square law force [9]. 
Consider spins 5*1 and 5*2 separated by a distance r that are at rest with respect to 
the preferred frame. Taking the Fourier transform of Eq. (8.21) with a; ^ 0, the 1/r 
potential between them takes the form 

T\/r2 1 

= {(9' + ■ 5^2 + (/ - l/a){S, ■ f){S, ■ r)) . (8.26) 

This reduces to the pure Goldstone result in the g' — > hmit, modulo a redefinition 
M \faM. This potential has an interesting feature which distinguishes it from 
both electromagnetism and ghost condensation. Consider a toroidal solenoid filled 
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with a ferromagnetic material.^ When current runs through the solenoid, the ferro- 
magnetic spins will align in the azimuthal direction, but as we will show in App. F, 

— * 

there will be no net magnetic field nor net longitudinal A field outside of the solenoid. 

— * 

However, there will be transverse A fields which can interact with a test spin, and 
because there is no magnetic field leakage from this geometry, this allows for a null 
test for the inverse-square law spin-dependent force. 

In order for the potential in Eq. (8.26) to be experimentally testable, it has to 
be at least roughly of gravitational strength, and if we assume spin sources with one 
aligned spin per nucleon, M/F has to be of order rrii^/Mpi ~ 10~^^. As shown in 
Ref. [9], M/F cannot be much larger than 10^^^ without the theory exiting its range 
of validity. When M/F is so small, it is easy to satisfy the condition in Eq. (8.25). 
In the previous section, we argued that the ratio of the vector coupling to the axial 
couphng can be made as small as mlMpi/ldn'^ ~ 10~^^ without fine tuning if we 
assume the quasi-parity invariance in Eq. (8.13). This tells us that 



< 10-14 =^ q^< 10-^5 (8.27) 



M/F 

in agreement with Eq. (8.25). 

This spin-spin potential is even more interesting at finite velocity. For sources 
moving at velocity v will respect to the preferred frame, the potential at fixed co- 
moving distance is given by the Fourier transform of Eq. (8.21) with the replacement 
uj ^ k ■ V. In the case f < c^, we can calculate the spin- spin potential as a power 
series in v/cg. For simplicity, we define 

f{A,B)^A-B-{A-f){B-f), (8.28) 

and the potential takes the form 

(8.29) 

- ^^^^5 {f{s,j2)mi)+2f{s,,v)ns2,v))+o{vyct). 

The unique angular dependence of this spin-spin potential is a smoking gun for gauged 
ghost condensation. As v gets larger than Cg, the potential should map onto the 
shadow/shockwave form of Ref. [9], though to see this behavior, one would need to 
keep track of k'^ terms in the dispersion relation for the Ai longitudinal mode. 



^Thanks to Blayne Heckel for suggesting this geometry. 
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Next, we consider energy loss due to the emission of Ai particles. In Ref. [9] we 
did not consider bounds on new light particles from astrophysics [57], because the 
scale of spontaneous Lorentz violation M was constrained to be much lower than 
typical astrophysical energies. Now that we can raise the scale of Lorentz violation 
well above even the electroweak scale, these bounds are relevant, but its easy to see 
that once we impose the experimental bound on the astrophysical bound is almost 
automatically satisfied. For a pseudo-scalar (p coupling to the axial current, stellar 
energy loss arguments bound the coefficient of the operator 

ij^aV (8.30) 

to be F > 10^ GeV [58]. In the non-rclativistic limit, Eq. (8.30) is identical to the 
canonically normalized coupling of the longitudinal mode of A to spin-density, so as 
long as M > 10 eV, the bound on is stricter than the bound on F directly. For the 

— * 

transverse modes of A, we can estimate the bound on M/F by going to canonical 
normalization and replacing the derivative with a typical astrophysical energy 
Castro ~ 100 keV: 

^4d^^r^^P-A^ => f9<10~''. (8.31) 

As long as M > 10~^ eV/g, the bound on fx enforces the astrophysical bound. This 
bound also applies to the vector coupling, requiring 

q^g < 10-^\ (8.32) 

Finally, for sources moving with respect to the preferred frame, there is both the 
possibility of Cerenkov radiation from gravitational couplings (sec Sec. 7) and from 
the coupling in Eq. (8.16).^° Because the scalar mode in gauged ghost condensation 
has a finite velocity, this Cerenkov effect is turned off when the velocity of the source v 
is smaller than Cg. When v/cg^ l,we expect the discussion of Ref. [9] to go through 
virtually unchanged, where wc found that within the range of validity of the effective 
theory, there were no observable effect from ether Cerenkov radiation for electrons. 
Ref. [59] studied the effect of Cerenkov radiation on neutrinos from SN1987A, where 
the bound was expressed as: 

Assuming that neutrino couplings are comparable to electron couplings, this does not 
place any additional bounds on couplings between A^ and axial currents. 

^°Note that there is no Cerenkov radiation from the couphng in Eq. (8.5), because it only involves 
couplings to the transverse modes of Ai, which travel at Cv, i.e. nearly the speed of light. 
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9 Conclusions 



We have analyzed spontaneous breaking of Lorentz symmetry by a vector condensate, 
a theory which gives rise to 3 extra modes in the gravitational sector. We showed that 
this theory can be viewed as a gauging of ghost condensation, the minimal theory 
of spontaneous breaking of Lorentz symmetry. The relation to ghost condensation 
allows us to understand the power counting of operators in the effective theory, and 
also to understand the limit where the gauge coupling g is small, where we recover 
ghost condensation. We also analyzed various constraints on this model, taking into 
account important nonlinear effects that are closely analogous to those found in ghost 
condensation [9]. 

We find that the limits on the scale M where Lorentz symmetry is spontaneously 
broken allow values much higher than the weak scale, suggesting that Lorentz viola- 
tion may be more closely connected with fundamental particle physics. Prom purely 
gravitational constraints, we find that M < Min(10^^ GeV, 10^^ GeV) from ob- 
servation of a steller-mass black hole and constraints on the PPN parameters. If we 
allow all possible couplings to the standard model and cut off gravity loop effects 
at the scale M, we find that we need M < 10^^ GeV from the absence of observed 
Lorentz violation. 

It is our hope that this work will bring Lorentz violation closer to the mainstream 
of particle physics, and that it will stimulate further investigations on the possible 
relation to the puzzles of gravity and cosmology. 



Alternative theories of gravity with additional vector fields pointing in some preferred 
direction have been considered in the literature [11, 12, 13, 14, 15, 16, 17, 18, 19, 21, 
22]. The value of the vector field is often imposed as a constraint through a Lagrange 
multiplier 



We will show that this corresponds to certain decoupling limit of our theory. To see 
this, we can add a term —2k'^/M'^ to the Lagrange multiplier. Integrating out k we 
obtain 



A Decoupling Limit and Compcirison with the Literature 



C D -Kiu^'u^ - 1). 



(A.l) 




(A.2) 
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which is exactly the leading term of our gauged ghost condensate Lagrangian: 



8 V M2 

making the appropriate mapping between u^^ and A/^. The constraint of the La- 
grange multiplier is obtained by sending M oo. Therefore, we sec that these 
theories correspond to the decoupling limit M — > cxo of gauged ghost condensation, 
with appropriate rescaling of other parameters, g~^,a, (3,'y, . . . while keeping 
g-^M^, aM\ (3M^, 7M2, . . . finite. 

In this appendix we will study this decoupling limit and see the extent to which 
they are healthy modifications of gravity. We start with the ungauged ghost conden- 
sate with the lagrangain 

C^jCeu + lM\X - 1)2 - lM\V\f, (A.4) 

where jCeh is the Einstein- Hilbert Lagrangian. For convenience we have defined — 
ccM^. We want to understand the physics of taking M — > 00, with M and Mpi held 
fixed. This is equivalent to imposing the constraint X — 1. 

The linearized Lagrangian for the scalar sector including a matter source p is 

£ = -M|i(V$)' + iM^($ - 7r)2 - lM\V\f + p*. (A.5) 
For large M, we impose the constraint $ = tt, which gives the effective Lagrangian 

= -M|,(V7r)2 - ^M^^^nf + pn. (A.6) 
From this we easily read off the dispersion relation for the scalar mode: 

Note that this has the wrong sign for > 0, which is just a manifestation of the 
Jeans instabihty in ghost condensation. But for < 0, the scalar has a healthy 
dispersion relation. 

Of course, the same results can be obtained from the dispersion relation for tt with 
finite M, Rewriting in terms of M: 

co^ = H -k^. (A.8 
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In the limit M — > oo only the first term survives. So the decoupling limit is equivalent 
to considering 

\k\^m=^ . (A.9) 

Therefore, we can define a healthy theory as long as m is large enough to be a UV 
scale and if we change the sign of relative to the ghost condensate. 

What is the Newtonian potential in this limit? Define the decoupling limit sound 
speed as 

In the \k\ <^m hmit, the propagator is 

At very late times (tc; — > 0), the theory does not have a Newtonian potential, so the 
decoupling limit of ghost condensation is not a viable modification of gravity at very 
late times. However, for times t short compared to r/vo, \k\ ^ u ^ vo\k\ 

M) ^ L^l. + o ( ^1 , (A.12) 

and we recover an ordinary Newtonian potential plus small corrections. If vq is 
sufficiently small, then the decoupling limit could be a viable modification of gravity 
with interesting late time behavior, and we will study this limit in future work. 

Now consider the gauged ghost condensate. To make the limit clearer, we define 
7r to have mass dimension —1, to have mass dimension +1, and rescale the fields 
and couplings relative to Eq. (4.3) by 

M ^ /MM 
/X \m) g^' 



- -l^iTj] -1^ (A-13) 



where /x is a mass scale held fixed when we take the limit M — > oo. The gauge 
transformations can now be written as 

5A,^d^X, Sn^-^, (A.14) 

A* 



and the scalar sector Lagrangian is 

1 



jC = -M|i(V$)^ + ^(VAof + lM\^ - TT + Ao/i^)^ - iM^V'nf. (A.15) 
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We are interested in the limit M — > oo with the other couphngs (including /j.) held 
fixed. The dispersion relation for the scalar mode is 

As before, the second term is absent in the limit M — > oo, and is negligible for 
\k\ <^m. The term has the right sign for 

or for < g < Qc- It is easy to see that in order to get a viable Newtonian 
potential at very late times {u — > 0), we must choose the first option. The static 
limit of the propagator is 

so if 5f < Qc this would have the wrong sign relative to the standard Newtonian 
potential. So as long as g > gc, the decoupling limit of gauged ghost condensation 
is a healthy modification of gravity at leading order. Interestingly, if < and 
g < gc, then there is once again an intermediate range of times for which there is an 
ordinary Newtonian potential. 

This analysis shows why, for example, Ref. [33] finds finite corrections to a^^^ in 
a model with the constraint in Eq. (A.l). In Sec. 4, we found 

which would diverge as oo. However in the decoupling limit where g^'^M'^ 

and aM^ are held fixed, a^^^ stays finite. We emphasize that while the decoupling 
limit is healthy, it obscures the power counting of gauged ghost condensation, which 
gives a systematic way to understand the relevance of operators in a Lorentz- violating 
setting. 



B Stress-Energy Tensor of the Gauged Ghost Condensate 

The goal of this appendix is to calculate the stress-energy tensor of the gauged ghost 
condensate. For this purpose we need to consider general variations of the metric 
since the stress-energy tensor is defined as the functional derivative of an action with 
respect to the metric components. On the other hand, in order to make the expression 
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independent of the heavy modes we have to integrate out the perturbation of Do(f> 
before taking the variation and, thus, we need to separate the time coordinate and 
space coordinates. For these reasons, we shall adopt the ADM decomposition of the 
metric: 

ds^ = -N^df + Qij (dx' + (3'dt) (dx^ + (3^dt) (B.l) 

and define by 

L>o0 = (1 + Ao)N, Di(t> = Ai. (B.2) 

We shall integrate out. 
Expanding X — 1 as 



\{X-l)^Ao-\ - ^] AA, + 0{Al) 



(B.3) 



the leading Lagrangian is 



1 
2 



A - ^ ( q'' 



^-^\aa- 



where g*-' = {q and q — det q. Integrating out the massive mode ^o, 



(B.5) 



we obtain 



F^^F'^" = -2q'^FuF^, + q^F^jF^i, 



(B.6) 



and 

where 
and 



1 

7^ 



d± { ^/q 



(3^Aj 
N 



kl 



+ 7^ U - 



AkAi 



+ -^^d,{N^q'^A,), 
(B.7) 

(B.8) 



Fm = 



1 

TV 



-{Ft^-P'Fk^) 



1 

2 



[3^ 

+ j^diAk, 



(B.9) 
(B.IO) 
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All components of the stress-energy tensor are given by taking the variation of 
the effective action w.r.t. N, and qij. The result is 



1 + - I g'' + 



Ar2 



»4.i;«4.i 



1 + 



N 



M2 



9 

1 - 



q''F,^F^, + f^A^Aj^dki^q'^F^i) 



1 + ^1 A5x(V^/^M 



N 



N 



kl 



(3^ 



rptj 



dtd^xNy/qjC 



^ { (^g''F±ikF±/ - Jg'™g'"Ffe/Fn^n) q'' - q%''F^kF±i + q'^q^'q^^FkmFin 



-q"'q''AkAi^dm{VQq"'''F^n) 
\1 



P'Ak 1 / 
N 2 



qkl_ 



AkAi \ d^{WD^ct>) + q'^'AkdiiWD^cf)) 



+q'''q^'AkAid4W^D,<P) - qViAkdi + Aidk){W^D,<P)} . 
The Einstein equation is 



where 



dt) ^ \dx^ ' 



(B.ll) 

(B.12) 
(B.13) 



Just for completeness, the equation of motion for Ai is 



•^3 AT 



1 



^dki^qq^'F^i) 
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+a 



N 



N ^ N J } 



(B.14) 



C Static Post-Newtonian Gravity 



The PPN parameters /9ppn and 7ppn are defined in the 1 / r-expansion of a spherically 
symmetric, static metric in the isotropic coordinate as 



2m 2/3ppNm2 3 3 

1 1 h 0{m /r ) 



27ppNm 2/ 2n 

1 H h 0[m /r ) 



(C.l) 



In this coordinate, /3ppn and 7ppn, respectively, measure the amount of non-linearity 
and the amount of space curvature produced by a mass. The values in General 
Relativity is Pppjq = jppN ~ 1- Experimental limits on them are 



|/3PPN - 1| < 10-^ I7PPN - 1| < 10-^ 

In this section we calculate /9ppn and 7ppn for the gauged ghost condensate. 
For the isotropic coordinate 

ds" = -N{rfde + B{rf {dr^ + r^dQj) , 

Einstein tensor is given by 



(C.2) 



(C.3) 



Cm = ^ 



^rr 



]_ (B^ 

G±_r — G±o — 0, 

1 



AB' 2B" 



rB 



B 



B' N' B'\ 2B' 2N' 

2 + — + + 

B N \B r B r N 



0. 



B^ B^_ (B^ 

'n'^ Vn^ Vb^ ^ ~ \ B . 



(C.4) 



We shall apply the formulae obtained in App. B to this metric by setting 

N = N{r), (3' = 0, Qijdx'dx^ = B{rf [dr^ + r^d^l) , Aidx' = A{r)dr. (C.5) 
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The result is 



9' 



1 /F, 



2\B 



1 + 



2^2 



a 



1,^,,2 A{dAY 



M-^T^e = 0, 



2B^ 



M-2rprr 



1 

1 



1 /F, 



2\B 



+ 



+ a 



r'^B'^ [2g^\B 



where 



N\l + 



2^2 



In order to calculate /3ppN and 7ppn, we expand all variables by 1/r, 



oo 

n=l ' 

oo rj 

1 + 



iV(r) 
5(r) 

n=l 

oo \ 

Mr) = E ^- 

n=0 ' 



(C.6) 



(C.7) 



(C.8) 



Accordingly, we obtain the 1/r-expansion of the Einstein equation. 

The (± r) -component of the Einstein equation is a{dAy — 0, which with the 
above 1/r-expansion implies that ^4 = 0. In this case non- vanishing components of 
the stress energy tensor are 

2 / 



9' 



ifN'Y ( A^ 



^2j^-2rprr ^ 



1 

'7 



2 \NB 



1 / N' 



r'^B-'^M-'^T^ 



J_ ( N' 



2 \NB^ 

2 



+ 



2^2 ; r^B^ \ 

^2 fr^BN'V 



'r^BN'' 
N 



2^5 



N 



(C.9) 
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The leading term in the (rr)-component of the Einstein equation says that 



(C.IO) 



With this relation, the leading term in G±± — B'^G"' = k^{T±± 



) is 



2g^Mp^, we obtain 




(C.ll) 



Hence, unless = 




(C.12) 



From (C.IO) and (C.12) we obtain 



/3ppN — 7ppn 



1. 



(C.13) 



Therefore, there is no constraint from the spherically symmetric, static post-Newtonian 
gravity. The leading correction to General Relativity appears in B2/NI. Actually, 
the leading term in the (_L_L)-component of the Einstein equation says that 



and this value of B2/NI has a deviation from the GR value 1/4. However, this 
deviation is in the post-post Newtonian order and beyond the current ability of gravity 
experiments. 



In Sec. 5 we have provided both analytical and numerical evidences showing that 
there should be no instability except for the perfectly plane symmetric case. We have 
also pointed out that, on the other hand, it is in principle possible to create a caustic 
plane if the perfectly planar symmetry is assumed. However, this situation is unlikely 
to happen in generic situations since the perfect plane symmetry is too ideal an 
assumption. Indeed, small fluctuations on top of the plane symmetric collapse should 
grow and the layer should be broken into pieces before a planar caustics occurs. After 
the fragmentation, a lower-dimensional caustics does not occur since, as discussed in 
subsection 5.3, for codimension 2 or higher it would cost infinite amount of work to 
compress a finite- volume, negative r^E,. [s = 1,2) region to an infinitesimal volume 
or thickness. 




(C.14) 



D Fragmentation of Plancir Caustics 
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In this appendix we perform a perturbative analysis supporting this picture of 
fragmentation. In particular, we analyze inhomogeneous linear perturbation on top 
of a perfectly planer caustic solution and show that inhomogeneous perturbation 
grows faster than the background caustics. 

In a fixed flat spacetime background ds^ — dt^ — dx^ — dz^, we consider a conflg- 
uration depending on two space coordinates x and z and the time t: 



Aidx' = Ax{t, X, z)dx + A^it, x, z)dz. 
The equation of motion is then 



(D.l) 



E, = A^ffE^ - d'Fij + ag^did'Aj, 
Ei = A- A^diAj, 



(D.2) 



where i and j run over x and z. We consider a planar scaling solution with inhomo- 
geneous linear perturbation of the form 



" -t 



l + f{t)e 



ikx 



Ax = 5(^)e^'^ 



(D.3) 



where we consider / and g as small perturbations. Since we are interested in the 
behavior near the caustic plane z — 0,we have truncated the Taylor expansion w.r.t. 
z at the lowest order. In the flrst order in / and g, the equation of motion becomes 



k'z' + 1 



f+tk 



2^ , 2^ , ..,..,2 / 



+ ikz h ik 

ty —t 



9_ 9 

-t ^ {-ty 

2^2 



0, 



(D.4) 



By eliminating g we obtain the third order equation for / as 



f 



-t 



(D.5) 



The general solution to this equation is 



k cos{kt) H sm{kt) 

t 



ksm{kt) -cos{kt) 



(D.6) 



where Ci, C2 and C3 are constants. Hence, the linear perturbation is unstable and 
grows faster than the background. 
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E A Simple Way to Calculate the Energy Loss in Cerenkov Radiation 



Consider a Lagrangian in the presence of an external source 

>CTotai = del)) + £5(0, S), (E.l) 

where Cs represent the source term, e.g. Cs — Sep. If the action without the source 
S(j, — J d^xC^ is invariant under the transformation — > 0' = + 5(f), there is an 
associated conserved Noether current 

with given by 

5C^ = d^V^". (E.3) 

The presence of the (fixed) external source "breaks" the invariance of the transfor- 
mation, so the current is no longer conserved: 

d,r^ = 5Cs. (E.4) 

Integrating over space, we obtain 

j (fx SCs = J d^x d^J^ = doQ^ + j rdSi. (E.5) 

The right hand side is just the total charge flowing into the (p field and out to infinity, 
so it must be compensated by the loss of the total charge of the source, —dQs/dt: 



dQs 
dt 



J d^xSCs- (E.6) 



For energy, the corresponding symmetry transformation is time translation, so 50 = (f). 
The energy loss for a source Cs = S(p is then simply given by 



dEs 
dt 



-jd^xS4). (E.7) 



In the linearized approximation, can be easily solved in terms of a given source 
in momentum space. It is then straightforward to compute the energy loss using 
Eq. (E.7). It applies to both the Cerenkov radiation when a source moving in a 
medium faster than the sound speed, and multipole radiations due to accelerations. 
As an example, we derive the energy loss from the Goldstone boson radiation in a 
binary system. The results are used in sec. 7 to derive bounds on the scale of the 
Lorentz violation from the binary pulsars. 
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We consider a system of binary pulsars where two pulsars of equal mass Mq move 
in a circular orbit separated by a distance 2ro. We ignore the velocity of the center 
of the system for this calculation. 

In the 5-function approximation, the source is given by 

p{x, t) = MoS^ {x - f{t)) + MoS^ {x + f{t)) , (E.8) 
f{t) — {rocoscuot, rosincuot, 0), (E.9) 

where the angular frequency ujq is related to the orbital velocity Vq by Vq = ujoVq. Its 
Fourier transform is 

p{u;,k)^2Mo e^'"(^'=-i)j2n(A;||ro)27r5(u;-2na;o), (E.IO) 



where = y/c^ + ky, t&nOk = ky/kx, and we have used the property of the Bessel 
functions, Jm{—x) = (— 1)™ J„i(a;)- A form factor f{k) can also be included to rep- 
resent the finite size effect. Applying the simple way of calculating the energy loss 
described above and using Eqs.(7.1)~(7.3), we obtain the time-averaged energy loss 
rate 



2 



J. , 2nuj, ^ 2 



(E.ll) 

Now we can consider two different limits. First, for vq = luoTo <C we can use 
the approximation 

1 / X \ ™ 

J^(x) ~ — r ~ ' small X. (E.12) 
m! V2/ 

The sum is dominated by the smallest \n\, n — ±1 which corresponds to quadrupole 
radiation. There is no form factor suppression (/(fc) — 1) because the inverse of the 
momentum is much bigger than the size of the system. The resulting energy loss by 
the quadrupole radiation is 

dEs ^ _4aM|MX ^ _2^^7mAPvl 

where in the last equality we have used the relation GnMq/tq = 4vq. 

In the opposite limit, Cg <C Vq, we need to use the asymptotic approximation of 
the Bessel functions. 



Jm{x) ^ \l — COS 



TTX 



X- (m + ^^^ , for X > 1. (E.14) 
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From Eq. (E.ll) the energy loss is given approximately by 

F A Null Test for a Spin-Dependent Force Law 

It is well known that in electromagnetism, the magnetic field outside of a divergence- 
free spin configuration vanishes. Similarly for (ungauged) ghost condensation in the 
non-relativistic limit, the coupling between spin and ghostone boson field vr takes the 
form [9] 

M 

Ant = -E^s- Vvr, (F.l) 

r 

and doing an integration by parts, we see that vr is not sourced by a divergence- 
free spin configuration. However, the potential in Eq. (8.26) is not of the form of 
Eq. (F.l). In particular, the transverse modes of Ai couple directly to s and not to 
its divergence. 

One interesting divergence-free spin configuration is a ring of spin. This could 
be established with a toroidal electromagnet, and because the magnetic field is zero 
outside of the torus, this allows a null test for the gauged ghost condensate spin 
potential. Consider the following setup, where a ring of radius R with net spin 5*1 is 
separated a distance r away from a point spin 5*2. 





It is straightforward to calculate the potential from Eq. (8.26): 

WS^fl^ (F.3) 

r ^ oirr r 

where S2y is the projection of 5*2 in the vertical direction. Note that this potential 
does not depend on a, because a only controls the coupling of the longitudinal mode 
of A. Though this potential goes as there is in principle no magnetic field 

leakage from the ring 5*1, so it should be easier to design a null experiment to look 
for this long-range spin-spin potential. 
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